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We study the current through molecular films as a function of orientation of the molecules in the 
film with respect to electrodes. It may change by more than an order of magnitude, depending on 
the angle between the axis of the molecules and the normal to the electrode. This is a consequence of 
a strong directional character of p-orbitals that determines the conductance through the molecules. 
We demonstrate this result on an exactly solvable model, and present the calculations for two 
different experimentally accessible molecular films sandwiched between gold electrodes. 

PACS numbers: 73.40. Gk, 73.61.Ph, 73.63.Rt 



The studies of organic molecules as possible electronic 
components and devices have grown out of the initial 
suggestion by Aviram and Ratner into a very active 
area of research, see e.g. One of the most impor- 

tant issues in this problem is the role of the molecule- 
surface interface. The geometry of the contacts [^|-^, 
the type of binding |H,0,D, and the molecule-electrode 
distance were all found to affect the conductance 

significantly. In this paper we perform a model study 
of the effect of mutual orientation of the molecules and 
the electrodes on conductance. We find that the ori- 
entational dependence of the current through molecular 
films is much stronger than reported previously in Ref. 

, and it should be taken into consideration in both in- 
terpreting experimental data and designing possible mo- 
letronic devices. 

The simple argument in favor of a substantial ori- 
entational dependence of the conductance is the large 
anisotropy of the molecule-electrode coupling. In most 
molecules that have been studied experimentally or the- 
oretically so far, the conduction is due to 7r-conjugated 
molecular orbitals (MOs). Such MOs extend over the 
whole molecule and facilitate the transport of electrons 
between the two electrodes. Since those MOs are made 
mostly of anisotropic p atomic orbitals, the overlap of the 
conducting MOs with electrode wave functions strongly 
depends on the angle between the main axis of the 
molecule and the surface normal. In general, we expect 
the overlap and the full conductance to be maximal when 
the lobes of the p-orbital of the end atom at the molecule 
are oriented perpendicular to the surface, and smaller 
otherwise, as dictated by the symmetry. An estimate 
of the effect can be made from the general properties of 
the p-orbitals. The overlap integrals of a p-orbital with 
orbitals of other types differ by a factor about 3 to 4 
for the two orientations jll]. Since the conductance is 
proportional to the square of the matrix element, which 
contains a product of two metal-molecule hopping inte- 
grals, the total conductance variation with overall geom- 



etry may therefore reach two orders of magnitude, and 
in special cases be even larger. This conclusion would be 
only valid when the electronic structure of the molecule 
and the molecule-electrode distance does not change sig- 
nificantly while its orientation is being changed. One 
can anticipate this to apply when the coupling to the 
electrode is weak, as might be the case in some experi- 
mental setups. The calculations (l2j suggest that there 
are attachements of particular molecules with very little 
resistance to the rotation with respect to metallic surface. 
The present effect (change of the geometry of the metal- 
molecule contact) should be distinguished from other ge- 
ometrical effects caused by miramolecular conformation 
changes. One possibility, discussed in Ref. is to ro- 
tate a segment of a 7r-conjugated molecular wire. That 
would break the conjugation and result in decreased con- 
ductance. Another would be to rotate an active part 
of the molecule with respect to the rest of it to effec- 
tively change the current path, which leads to significant 
changes in resistance of the effective circuit [Q. 

In order to illustrate the geometric effect on current 
we shall first consider a simple model shown in Fig. |l| 
(top panel). A molecular film is sandwiched between 
two three-dimensional electrodes, which are assumed to 
have a simple cubic structure with one s-orbital per site 
and the onsite energy es, which we take as an energy 
origin. We assume, without loss of generality, that the 
band in the leads is half-filled, so that the Fermi level 
Ep = 63 = 0. The hopping integral between the neigh- 
bors is —ts- The "molecules" have only two p tt— type 
atomic orbitals with the onsite energies Cp and hop- 
ping between those tt— orbitals is —t-^. Those orbitals 
model TT-conjugated bands of organic conductors. The 
molecules are allowed to rotate about the bottom "atom" 
and thereby change the angle 9 between the main axis of 
the molecules and the surface normal. As 9 changes, the 
p-orbitals remain oriented perpendicular to the molec- 
ular "backbone". Because of that, the matrix element 
between a molecular p-orbital and an electrode s-orbital 
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FIG. 1. Top panel: the molecular film between two leads 
with simple cubic structure and one s-orbital per site. Only 
one molecule of the film is shown explicitly. Due to the sym- 
metry of coupling to the electrodes, the current is zero for 
the normal orientation (d = 0) and maximal for the parallel 
orientation (Q = it/2) of the molecules. Bottom panel: cur- 
rent per surface unit cell versus bias voltage for the tilt angles 
9 = 7r/2, 7r/4, and n/8, with the hopping parameters t^, t^^, 
and tsp in arbitrary units (to). Bias is in units of io/|e| , and 
the parameter Jo — 10~'^io^. 



varies as tspSinO, where t^p is a constant {sp— hopping 
matrix element between the molecule and the electrode). 
The molecule-electrode bond lengths do not change with 
6, and the surfaces of the two electrodes always remain 
parallel. Each molecule is assumed to have one electron 
per atom. The distance between the lowest unoccupied 
MO (LUMO) and highest occupied MO (HOMO), or the 
"HOMO-LUMO" gap, is Ei^umo - Enouo = 2t,. The 
middle of the gap is aligned with the chemical potential 
of the leads at zero bias voltage V = (meaning that 
ep = Cs = 0). 

The current through the film at low temperatures is 
given by a standard expression [|l5[|l6| ] 
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where the summation goes over the surface Brillouin zone 
of the lead for open channels, and the transmission co- 



efficient is found from solution of the scattering problem 

t{E, fc|| ) = V^^I^t.t^pe^('=^^+'=^«) sin^ e, (2) 
where 

42 ^^.tk^L +£*'=-«) sin^ 6* 

(3) 
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and we assume Ti = 1. Then, the band velocities are 
VzL(R) = '^tsShikzL(R), with kzL(R) the z-components of 
the momenta in the leads, which are found from 
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It is instructive to analyze the transmission in the limit- 
ing case of weak coupling to the electrodes, t^p/t^ <C 1. 
The exact transmission probability, Eq. (H), acquires 
the resonant form for zero bias voltage, V — Q, when 
kzL = kzR-. 



with 



£^1(2) = +(-)ix 



A = A (^S, fc|| ; F j = [tip/ts) cos k, sin" 0, 
r = r (S, fcji ; V) = {t%/ts) sin kz sin^ 9, 
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where A, F ^ i?i — i?2 = ^t-r^ is the shift of the resonant 
levels, and their width, respectively. For non-zero bias 
the transmission probability is characterized by different 
line widths for hopping to the left (right) lead F^, (F^) , 
and Eq. (||) will have a general form X]r=i,2 (g-g^)^+rV4 
with F = Fl -t- Ffl. If we were to neglect the dependence 
of F on energy and k\\ , we would obtain 
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The low bias case corresponds to sub-resonant tunneling 
inside the HOMO-LUMO gap, whereas the case of large 
bias involves the tunneling through HOMO and LUMO. 
There is a crossover in the angular dependence of the cur- 
rent through the molecule from sin* 9 to sin^ 9 with the 
bias. The current saturates at large bias (with further 
changes due to the density of states in the leads) and its 
angular dependence becomes less steep. This behavior is 
confirmed by the results of exact calculations shown in 
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ric position of the molecules with respect to the surface. 
Those molecules attach strongly to the gold substrate 
by thiol, -S-, end groups, which form a chemical bond 
with Au 18 1 . The strong bonding might well mean 
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FIG. 2. Top panel: schematic of benzene-(l,4)-dithiolate 
molecules positioned between two gold electrodes. Also shown 
is the LUMO of the free molecule found from the density 
functional theory. Bottom panel: conductance as a function 
of the tilt angle 6. 



Fig. It shows a dependence of current on the angle, 
which is stronger at small biases. The maximum of the 
current corresponds to 9 = 7r/2, i.e. to the horizontal 
orientation of the molecule in the film (backbone of the 
molecule parallel to the surface of the electrodes). Ex- 
act results also show the non-monotonous behavior of the 
current with the bias, when it is comparable to HOMO- 
LUMO gap, which comes from the complex bias depen- 
dence of k^^ fj^ through the dispersion relations (^. 

The formula (||) is, obviously, a particular form of the 
Breit-Wigner formula and it should apply whenever the 
width of the resonances is much smaller than the separa- 
tion between them. Thus, the results, Eq. (^, are more 
general and should qualitatively apply in more complex 
molecules too, if the overlap between resonances is not 
substantial, although it is likely to be violated at high 
bias voltages. 

We illustrate the angular dependence on two types 
of molecules, benzene- (l,4)-dithiolate (-S-C6H4-S-) and 
a, a'-xylyl-dithiolate (-S-CH2-C6H4-CH2-S-) placed be- 
tween two gold electrodes. Both of those systems 
were studied experimentally ||l^,|l8| and theoretically 
10|jig|,|2Cl[| , but mostly assuming a particular symmet- 



a substantial change in the electronic structure of the 
molecule, because of a possible charge flow to/from the 
molecule bonded to a metal (electron reservoir) due to 
the difference in their ionization potentials. Therefore, 
the following calculation is rather the model study to il- 
lustrate the orientation effect, which may not be as strong 
in a particular setup compared to our calcualted 

effect. It is clear, however, that the effect exists and may 
be large in situations when the coupling to electrodes 
does not strongly perturb the molecule itself. 

In order to compute the conductance, we follow the 
same general procedure of Rcf . . The gold electrodes 
have been described by a tight binding model with nine 
s— , p— , and d— orbitals per each Au atom with the pa- 
rameters from Rcf. The equilibrium molecular ge- 
ometry is found by the total-energy density functional 
minimization ||2^. The tight-binding parameters for the 
molecules and molecule-lead interfaces have been taken 
from the solid-state table of elements . 

Consider the benzene-(l,4)-dithiolate film. We first 
place the molecules perpendicular to the gold electrodes, 
as shown at the top of Fig. ||. Also shown is one of the 
conducting orbitals, the LUMO. For this orientation the 
conductance should be small because the p-orbitals, that 
constitute the LUMO (and the HOMO, which also con- 
ducts), are parallel to the surface and the overlap is the 
smallest. Next, we allow the molecules to rotate about 
one of the end sulfur atoms while keeping the two surfaces 
parallel and both sulfur- gold bond lengths constant. The 
molecules are rotated in such a way that the lines through 
pairs of hydrogen atoms, that are symmetric with respect 
to molecules' main axis, remain parallel to the gold sur- 
faces. As the angle 9 between the molecules and the 
surface normal increases, so does the overlap of the sul- 
fur p-orbitals with the gold orbitals. Thus, we expect the 
conductance to grow with 9. 

Our numerical results are shown in Fig. ^ We have 
studied two different surfaces, Au(lll) and Au(lOO), and 
two binding schemes, the "on-top" and the "hollow" po- 
sitions. In the first case, the sulfur atoms are positioned 
directly on top of the gold atoms at a distance of 2.39 A 
|l2| . In the hollow case, the sulfurs are in symmetric posi- 
tions with respect to three gold atoms of the Au(lll) sur- 
face, and with respect to four gold atoms of the Au(lOO) 
surface. In all four cases, we found the conductance to be 
a steadily increasing function of 6. The ratio G'(-|)/G(0) 
is the highest (=26) in the Au(100)-hollow case and low- 
est (=17) in the Au(lll)-hollow case. The shapes of all 
curves are qualitatively similar, while the hollow posi- 
tion is 2 to 3 times more conductive than the top posi- 
tion. Thus we observe a large orientational effect, which 
magnitude is rather insensitive to the electronic structure 
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FIG. 3. Top panel: schematic of a, a'-xylyl-dithiolate 
molecules placed between two gold electrodes. Also shown 
is the LUMO of the free molecules computed within the den- 
sity functional theory. Bottom panel: conductance of the film 
absorbed on the Au(lll) surfaces. 

and geometry of the electrodes. This finding corroborates 
our main argument that the orientation of molecular p- 
orbitals with respect to the electrodes is an important 
factor in the problem. 

Next, we consider another molecule, a, a'-xylyl- 
dithiolate, see Fig. ^. The 9 = position is chosen to be 
the one with the S-C bond perpendicular to the gold sur- 
face. Then, again, the angle between the S-C bond and 
the surface normal is systematically changed in both di- 
rections, while keeping the gold-sulfur bond lengths con- 
stant and the two electrode surfaces parallel. Since this 
molecule is asymmetric with respect to the S-C bond, 
the allowed interval of the tilt angle is also asymmet- 
ric, —OAn < 9 < 0.2tt (positive angle corresponds to the 
clockwise rotation of the molecules) . We have studied the 
Au(lll) surface and the two binding positions of sulfur, 
the on-top and the hollow. The computed conductance 
is shown at the bottom of Fig. |^. It is maximal at small 
9 and falls off upon tilting in either direction. Such a be- 
havior is readily understood if the spatial structure of the 
conducting MOs is taken into account. In Fig. || one of 
those orbitals, the LUMO, is shown. The hybridization of 
the sulfur p-orbitals with the rest of the molecule has led 
to the wave function having a preferential direction from 
the sulfur atoms towards the ring. Still, the wave func- 



tion retains its p-character so that a strong anisotropy of 
the LUMO-gold overlap is expected. Note that, accord- 
ing to Fig. I the precise direction of the wave function is 
in slight misalignment with the S-C bond, which is cho- 
sen as 9 = 0. The maximum overlap, and consequently 
maximum conductance, is therefore expected for a small 
negative 9, precisely the behavior we obtained from the 
full numerical treatment of the scattering problem. This 
detail indicates the adequacy of the present tight-binding 
parameterization. 

The present results show that the orientational depen- 
dence "molecule-electrode" of the conductance is a large 
effect if the conducting orbitals are extended 7r-orbitals, 
and the contact with electrode does not strongly perturb 
the molecule itself. The effect derives from the strong 
directional character of 7r-orbitals. We have also studied 
several other systems of different complexity, both ana- 
lytically and numerically, and obtained the results con- 
sistent with this simple picture. It is usually assumed 
that the geometry in molecular transport experiments is 
fixed, but there may be circumstances when the geome- 
try can be changed. Firstly, for some molecule-electrode 
pairs there may be more than one possible chemisorp- 
tion mode. For instance, for thiolates on gold it was 
found in Ref. that the hollow and on-top binding 
configurations are close in energy but have very differ- 
ent surface-S-C bond angles, ^ 180° and ^ 105°, re- 
spectively. Consequently, the molecules may change the 
configuration with temperature and/or under other ex- 
ternal factors. According to our calculations, this may 
lead to a significant change of conductance. Secondly, 
the tilt angle in self-assembled monolayers or Langmuir- 
Blodgett films, used in moletronics studies, is never ex- 
actly zero [ p3[ | and may depend on the temperature, pres- 
sure, and/or other parameters of preparation. Thirdly, 
the angle may change during the measurement process. 
For instance, in STM probes of single molecules the 
position of the probing tip is not fixed with respect to 
the molecule and may record a distribution of conduc- 
tance if the latter is strongly angle-dependent. Finally, 
a finite temperature always results in geometry fluctua- 
tions that in turn lead to conductance fluctuations. Thus 
the present mechanism may be one of the sources of the 
temperature dependence of the conductance. 
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